Spinon heat transport and spin-phonon interaction in the antiferromagnetic spin-1/2 
Heisenberg chain cuprates Sr2Cu03 and SrCu02 
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We have investigated tiie tliermal conductivity K,nag of liigli-purity single crystals of the spin chain 
compound Sr2Cu03 which is considered an excellent realization of the one- dimensional spin-1/2 
antiferromagnetic Heisenberg model. We find that the spinon heat conductivity Kmag is strongly 
enhanced as compared to previous results obtained on samples with lower chemical purity. The 



analysis of 



allows to compute the spinon mean free path Zmag as a function of temperature. 



At low-temperature we find /mag ~ 0.5 pm, corresponding to more than 1200 chain unit cells. Upon 
increasing the temperature, the mean free path decreases strongly and approaches an exponential 
decay ~ i exp T^ /T which is characteristic for umklapp processes with the energy scale ks Tu . Based 
on Matthiesen's rule we decompose Zniag into a temperature-independent spinon-defect scattering 
length Iq and a temperature dependent spinon- phonon scattering length lsp{T). By comparing 
^mag(T) of Sr2Cu03 with that of SrCu02, we show that the spin-phonon interaction, as expressed 
by Isp is practically the same in both systems. The comparison of the empirically derived Igp with 
model calculations for the spin-phonon interaction of the one-dimensional spin-1/2 XY model yields 
reasonable agreement with the experimental data. 

PACS numbers: 75.40.Gb, 66.70.-f, 68.65.-k, 75.10.Pq 



I. INTRODUCTION 



The physics of low- dimensional quantum magnets has 
recently attracted considerable attention by experimen- 
tally and theoretically working scientists because intrigu- 
ing properties are found. Such systems exhibit a variety 
of unusual ground states and exotic elementary excita- 
tions which, in some cases are well accessible by theoret- 
ical treatments. An important class of materials which 
host such quantum magnets is formed by copper-oxides 
(cuprates) which feature Cu^^-ions which, through their 
3(P configuration, generate S = 1/2 sites. The type and 
strength of the interaction between these spins depends 
crucially on the structure of the material. For example, 
amongst the cuprates there are model systems which real- 
ize S = 1/2 arrangements with strong antiferromagnetic 
(AFM) Heisenberg-type interaction {J/ks ^ 2000 K) 
in the form of square lattices, two-leg spin ladders and 
chains About 10 years ago it was discovered, that in 
all these different systems their magnetic excitations give 
rise to a highly anisotropic thermal conductivity tensor of 
the respective materials, with an unexpectedly large mag- 
netic contribution along the directions of large AFM ex- 
change. This finding opened up a new route to investigat- 
ing the generation, scattering and dissipation of magnetic 
quasiparticles (complementary to neutron scattering and 
magnetic resonance experiments) through analyzing the 
magnetic thermal conductivity Kmag^''— Among these 
findings, the results for compounds which realize the 



one-dimensional S = 1/2 AFM Heisenberg model (ID- 
AFM-HM) are particularly interesting because funda- 
mental conservation laws predict ballistic heat transport 
in these systemsJii^ This means that in the ID-AFM- 
HM model, an infinite magnetic heat conductivity is ex- 
pected. Despite this rigorous prediction, in any real sys- 
tem, the transport is dissipative, due to extrinsic scatter- 
ing mechanisms. Nevertheless, an unprecedentedly large 
Kmag has been observed in ultra-pure samples of the zig- 
zag chain compound SrCu02 J^^^ More specifically, upon 
enhancing the chemical purity of the compound, Kmag be- 
comes i) increasingly larger, and ii) even for the highest 
purity level it can be fully described (in the framework 
of a simple kinetic model) by considering spinons scat- 
tering off impurities and phonons only, where the impu- 
rity scattering fully accounts for purity dependence of 
Kmag and the phonon scattering prevails at elevated tem- 
peratures. This is indeed consistent with the predicted 
ballistic transport since no further scattering (i.e. spinon 
spinon scattering) mechanism needs to be invoked Fur- 
thermore, the analysis of the data yields very clean data 
for the spinon-phonon scattering, for which a full theo- 
retical description is still lacking. 

In the zig-zag chain compound SrCu02, two 5=1/2 
chains with large AFM exchange are tied together by a 
significant but frustrated intrachain exchange. In this 
paper we extend the previous findings for SrCu02 to the 
sm^Ze-chain material Sr2Cu03, where such complications 
are absent. We find that in high purity samples of this 
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compound Kmag is strongly enhanced as compared to pre- 
vious results^i^ for lower purity. Upon comparing the 
dependence of Kmag on the nominal purity level we ob- 
serve that Kmag of Sr2Cu03 depends in a similar fashion 
on the purity of the material like SrCu02. More specifi- 
cally, within a simple kinetic modetS, the spinon mean free 
path /mag can be decomposed into terms which describe 
spinon-defect and spinon-phonon scattering as is the case 
for the zig-zag chain compound. We show that using an 
empirical formula^ for describing the spinon-phonon scat- 
tering that the strength of this mechanism is practically 
indistinguishable for both materials. Going beyond this 
empirical approach we model the spinon-phonon scatter- 
ing by employing results for the spin-phonon interaction 
of the XY-mode\ which further underpins these findings. 
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II. MATERIALS DETAILS 
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The main building blocks of Sr2Cu03 are corner shar- 
ing chains formed by CuOa plaquettes<2i The chains 
are parallel to the crystallographic &-axis in Sr2Cu03. 
The intrachain exchange interaction between neighboring 
Cu^+ sites is mediated by 180° superexchange through 
the oxygen and with J/ks « 2150 - 3000 K^i^ is 
among the largest among known Heisenberg spin chain 
materials. The Cu-O-Cu chains are separated by Sr 
atoms, leading to an extremely small interchain inter- 
action of Jperp/ks ~ 0.02 K. Furthermore, muon spin 
rotation and neutron scattering measurements^iiS^ show 
that this material orders three-dimensionally only below 
the Neel temperature of « 5.4 K. These properties 
make Sr2Cu03 an excellent realization of a 5 = 1/2 
Heisenberg chain for temperatures T > 5.4 K. While 
Sr2Cu03 consists of isolated chains, the main structural 
element in SrCu02 is formed by Cu02 zig-zag ribbons, 
which run along the crystallographic c-axis. Each ribbon 
can be viewed as made of two parallel chains of corner- 
sharing Cu02 plaquettes, where the straight Cu-O-Cu 
bonds between corner-sharing plaquettes of each double- 
chain structure result in a very large antiferromagnetic 
intrachain exchange coupling J/fes « 2100 — 2600 K of 
the 5 = 1/2 spins at the Cu^+sitesj^i^^ The interchain 
coupling J between Cu^"*"-sites of two edge-sharing pla- 
quettes is much weaker (J /JwO.l — 0.2) 1^*21 Frustra- 
tion of this exchange interaction and presumably quan- 
tum fluctuations prevent three-dimensional long range 
magnetic order of the system at T > Tjv « 1.5 — 2K w 
lO^^J/fcs K3^i^ Hence, at significantly higher T the 
two chains within one double chain structure are usually 
regarded as magnetically independent. In fact, low-T 
(12 K) inelastic neutron scattering spectra of the mag- 
netic excitations can be very well described within the 
5 = 1/2 Heisenberg antiferromagnetic chain model<^ 



Figure 1: Thermal conductivity of Sr2Cu03 parallel to the 
spin chains along Hb for various purities. The dashed lines 
represent results from Sologubenko et. al. with 2N purity, re- 
produced from Ref.0. The dash-dotted line has been obtained 
by Kawamata et. al. for 3N purity and is reproduced from 
Ref . [20I. Inset: Thermal conductivity of Sr2Cu03 perpendic- 
ular to the spin chains along Kc for 2N (also reproduced from 
Ref. 0) and 4N purity. The solid lines are fits to the Callaway 
model. 



III. EXPERIMENTAL DETAILS 



Large single crystals of pure Sr2Cu03 were grown by 
the traveling solvent floating zone methodi^ The feed 
rods were prepared using the primary chemicals CuO 
and SrCOs with 4N (99.99%) purity. The crystals react 
rather rapidly with water. A hydroxide layer is formed 
at the surface of the crystals, after a short exposure to 
Therefore the crystals were annealed at 900°C in 
an oxygen atmosphere for three days before further treat- 
ment. This reverts the hydroxide layer back to Sr2Cu03, 
although in a polycrystalline form. 

The crystallinity and stoichiometry of all crystals were 
checked under polarized light and by energy-dispersive 
X-ray spectroscopy, respectively. For the transport mea- 
surements rectangular samples with typical dimensions 
of (2-0. 5-0. 5) mm'^ were cut from the crystals for each 
doping level with an abrasive slurry wire saw. Four- 
probe measurements of the thermal conductivity k were 
performed^^ in the 7-300 K range with the thermal cur- 
rent along the b and c-axes (k;, and Kq respectively). 
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IV. EXPERIMENTAL RESULTS 



Fig. [T] presents our findings for the lieat conductivity 
of Sr2Cu03, measured with the heat current parallel to 
the b and c axes, i.e. parallel and perpendicular to the 
chains in the material. We focus first on the temperature 
dependence of the thermal conductivity perpendicular to 
the spin chain, Kc, which is shown in the inset of the fig- 
ure. Along this direction, the heat conductivity of this 
electrically insulating material is purely phononic: As a 
function of temperature, it shows a characteristic peak 
at T = 22 K. and then strongly decreases upon further 
rising the temperature. The height of the peak sensi- 
tively depends on the density of impurities in the sys- 
tem, which generate phonon-defect scattering. This can 
be well inferred by comparing our data for a 4N-purity 
material with that of 2N (i.e., 99%) purity, taken from 
Ref. 1^. For this lower-purity sample the overall magni- 
tude of Kc is strongly reduced as is expected for typical 
phonon heat conductorsjSi In fact, the data for both pu- 
rities can be well described in the framework of a model 
by Callaway^, where the difference between both curves 
is largely captured by different point defect scattering 
strength (see Appendix). 

The thermal conductivity parallel to the chain, k^, is 
shown in the main panel of Fig. [TJ Kf, exhibts a peak 
at the same position as observed for Hc- However, the 
peak is much broader and the overall magnitude of K{, is 
significantly larger than that of the purely phononic Kc, 
which is the signature of a substantial magnetic contri- 
bution in this material, i.e., the heat current parallel to 
the spin chain is carried not only by phonons but also 
by spinonsi^ For our 4N purity sample, the anisotropy 
between Hb and Kc is roughly constant above 100 K and 
approximately Hb/Kc ~ 3.5. This is significantly larger 
than the previously reported^ anisotropy for 2N-purity 
Sr2Cu03 and provides clear evidence that the enhanced 
purity leads to a relative enhancement of the spinon con- 
tribution to the overall heat conductivity. The purity 
dependence of Kf, can directly be read off from the figure 
where we compare our findings with experimental data 
for 3N (99.9% purity) and 2N puritysamples of SraCuOs, 
reproduced from Ref. [^^ and Ref. H, respectively. From 
low to intermediate temperatures [7K — 150 K], the heat 
conductivity is strongly enhanced upon increasing the 
sample purity. At higher temperatures this pmrity de- 
pendence becomes weaker since the curves approach each 
other. Apparently it is possible to separate the thermal 
conductivity into two distinct regimes, where different 
scattering processes dominate. At low-T, the extreme 
sensitivity to impurities suggests, that spinon scattering 
off defects is dominating. At high-T, a further extrinsic 
scattering mechanism - spinon-phonon scattering^ - be- 
comes dominating as a consequence of increasing phonon 
population, which leads to the very similar Kf,^2N and 
KbAN for T > 200 K. 
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Figure 2: Estimated magnetic thermal conductivity of 
Sr2Cu03 for 4N (circles, squares), 3N (dash-dotted line)^, 
and 2N (dashed line)^ purity. The shaded area illustrates 
the uncertainty from the estimation of the phononic back- 
ground. The 4N results shown in squares instead of circles 
have a large uncertainty. The thick solid line is a fit for 
T > 80K with Kmag.flti ~ The thin solid line is a 

fit with «:niag,flt2 ~ exp {T*/T). 



V. DATA ANALYSIS AND DISCUSSION 
A. The spinon heat conductivity of Sr2Cu03 

The thermal conductivity parallel to the chain is com- 
posed of a magnetic and phononic contribution Kb = 
'^mag + 'i6,ph- In Order to extract the heat conductivity of 
the spin chain the phononic background is approximated 
as Kb^ph ~ Kc- This simple assumption is reasonable, 
since the purely phononic anisotropy between and Kc 
is small.i^ Therefore Kb^ph is not expected to be much 
different. The thus obtained spinon heat conductivity 
^mag = Kb — Kc for the 4N sample, as well as the results 
from literature for 3N and 2N, for which Kmag was ex- 
tracted in a similar way^ii^ii^ are shown in Fig. [2] Below 
T < 40K. i.e., in the vicinity of the peak of errors 
become large and the data in this range are neglected for 
further analysis. For higher T, the possible uncertainty of 
'*mag is around ±15%, which accounts for the individual 
errors of Kb and Kc- 

Starting from low-T, K,nag of the 4N sample increases 
almost linearly towards a peak. The uncertainty is 
quite large in the temperature regime, and hence an 
increase with a higher powe r^^i^^ of T in this regime 
(T < 40 K) cannot be excluded. However, in the sim- 
plest case of a temperature independent spinon-defect 
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scattering rate, one expects Kmag to be directly propor- 
tional to the thermal Drude weight I?tii, which also in- 
creases linearly with temperature at low-T up to T 
O.lSJ/fcs - 300 K,iii^-ii The peak is quite pronounced 
and found at ^48 K with a maximum value of about 
ISOWm^^K^^. The peak is followed by a strong de- 
crease for higher T. Such a temperature dependence can 
not be accounted for from the T-dependence of the ther- 
mal Drude weight, since it is expected to decrease only 
at very high T > O.GJ/ks ^ 1200 KM'^^ Instead, the 
decrease is consistent with our earlier notion of dominant 
spinon-phonon scattering at high-T. Kmag of the 3N and 
2N samples is qualitatively very similar to that of the 4N 
sample but is increasingly suppressed with growing impu- 
rity level, accompanied with a shift of the peak-position 



of Kn 



towards higher temperatures. Furthermore, the 



curves approach each other with increasing temperature, 
as the Kb data. 

Chernyshev and Rozhkov have proposed a model which 



describes 



^mag 



of Sr2Cu03 at the 2N purity level very 



^g^ ^36,37 However, a similarly convincing description of 
our data for the 4N sample is not possible. The same 
holds for the double chain material SrCu02- Spin- 
phonon drag has been suggested as a possible explana- 
tion for the failure of the model 4^ Recently there has 
been substantial progress in the formal theoretical treat- 
ment of this phenomenon)^ However, specific model cal- 
culations have not yet been performed for the materials 
under scrutiny here. It thus remains unclear whether 
spin-phonon drag plays a significant role in our experi- 
mental data, and we analyse Kmag without taking into 
account a possible contribution due to this effect. 

The temperature dependence of K^ag at high temper- 
ature T > 80 K up to room temperature can equally 



well be described by either k 



mag 



const, or 



'^mag cx exp(r*/T), with T* a characteristic energy scale. 
Fig. [2] shows the corresponding fits. While the physi- 
cal meaning of the former functional form remains elu- 
sive, one expects the exponential one for spinous scatter- 
ing off phonons from general considerations for Umklapp 
processes 4^ 



B. The spinon mean free path 



We proceed by calculating the mean free path of the 
spinons, imag, from the experimental Kmag byii^iiiiii^ii^ 



3^ 



7riV,,fc|T 



(1) 



where Ns is the number of spin chains per unit area. 
The thus extracted mean free paths are shown in Fig. [3l 
ImngiT) is very large at low temperature 0.5 jim) and 
decreases strongly with increasing temperature, consis- 
tent with the above already inferred increasing impor- 
tance of spinon-phonon scattering. While the high tem- 
perature data (r > 100 K) reflects well the exponential 

suppression of Kmag, i-e. ^mai 
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Figure 3: Derived mean free paths of the spinon excitations 
in Sr2Cu03 for 4N (circles, squares) purity. The solid black 
line is a fit to the data as explained in the text. The squares 
in the 4N data have a large uncertainty and are disregarded 
for the fit. A comparison is made with the fits for 2N (dashed 
line)*, and 3N (dash-dotted line)'^" purities. The dotted hne 
is a fit done to the 4N results by lap (eq. [2} only. The shaded 
area illustrates the uncertainty from the estimation of the 
phononic background. The inset shows the same results in a 
double logarithmic scale. 



a dotted line in Fig. [3]), the low temperature data clearly 
deviate from this functional form, which indicates that 
spinon-defect scattering becomes important. In order to 
test this notion, we compare in Fig. [3] these /mag data 
for the 4N sample with fits to the results of 2N and 3N 
purity samples as given in Refs. I8ll20l As expected, the 
enhanced impurity density in these samples causes a cor- 
responding reduction of Zmag at low temperature. Note 
that at high temperature the curves cling to that of the 
4N sample, which corresponds to the natural expectation 
of an identical spinon-phonon scattering strength in all 
samples. 

In order to capture this behavior in our further anal- 
ysis, we apply Matthiessen's rule for the scattering pro- 
cesses of the spinons l~ag (^) = ^a^+^sp (^)i where Iq de- 
notes the T-independent spinon-defect scattering, while 
Igp (T) takes the T-dependent spinon-phonon scattering 
into account. According to our empirical finding of an 
exponential decay of Kmag(r) at high temperature and 
consistent with previous findings^^^ we estimate Isp (T) 
as 



Y cxpT*/T, (shown as 



"sp 



exp(r*/r) 



(2) 
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Figure 4: Magnetic mean free paths of SrCu02 (pentagon 
shape) and Sr2Cu03 (circles) for 4N purity. The shaded area 
illustrates the uncertainty from the estimation of the phononic 
background. The lines are fits to the mean free paths. 



with fit parameters T* and Ag. As can be seen in Fig. [3] 
such a empirically derived functional form for the spinon- 
phonon scattering allows an excellent fit of the data with 
Iq = 0.54±0.05 ]jm (corresponding to approximately 1367 
lattice spacings), T* = 210 ± 11 K, and = (6.1 ± 5) x 
lOWiR-i. 



C. Comparison of SrCu02 and Sr2Cu03 

The drastic enhancement of Kmag and ^mag upon in- 
creasing the purity level of Sr2Cu03 provides strong ev- 
idence that the spinon heat transport in the S = 1/2 
AFM Heisenberg model as realized in this material is 
only limited by the external scattering off defects and 
phonons. Our findings here thus further corroborate pre- 
vious experimental evidencoi^ for the ballistic nature of 
heat transport in the S = 1/2 AFM Heisenberg model, 
obtained for the zig-zag chain compound SrCu02. Since 
the individual chains in SrCu02 and in Sr2Cu03 consists 
of the same structural elements (Cu02 plaquettes), it is 
instructive to directly compare the findings for spinon- 
phonon scattering obtained for both compounds. For 
this purpose we show in Fig. |4] the spinon mean free 
path /mag(r) of high-purity (4N) SrCu02 and SrsCuOa 
in units of the chain unit cells. In this representation, 
the mean free path is free of any geometrical particulari- 
ties and can directly be related to distances of spin sites 
within a single chain. As can be seen in the figure, the 
mean free paths for both compounds are virtually iden- 



tical for temperatures T > 150 K (the relative difference 
is less than 5%). In the whole temperature range, both 
curves can be well fitted with expression [5] for the spinon- 
phonon scattering, with the same T*, A^ as determined 
afore and with and Iq = 1.56 ± 0.16 jim (correspond- 
ing to approximately 3984 lattice spacings) for SrCu02. 
In fact, the same parameters for Igp (T) could be used 
for both compounds. This demonstrates that the spinon 
phonon interaction is the same in both compounds, de- 
spite the difference of their Cu02 chain structures. Since 
the chains in SrCu02 and Sr2Cu03 are composed of the 
same copper-oxygen plaquettes, one has to conclude that 
only phonon modes which modulate the Cu-O-Cu bonds 
along two corner-sharing plaquettes lead to a significant 
scattering of spinous. 

At low temperature, the mean free path of SrCu02 is 
by a factor of three larger than that of Sr2Cu03, in spite 
of the same nominal purity. The higher defect density of 
Sr2Cu03 may be caused by the relatively lower chemical 
stability. It reacts quite rapidly with water and decays, if 
exposed to air for a few hours. Additionally, with regard 
to unavoidable intrinsic crystal defects, the double chain 
is a much more stable structure, due to its layout of the 
Cu02-plaquettes. 

Apart from these minor and plausible differences with 
regard to the spinon heat conduction, we would like to 
point out a surprising dissimilarity which concerns the 
phonon heat conductivity Kph. As we have seen in Fig.[Tl 
Kph of Sr2Cu03 increases strongly at increasing the pu- 
rity, as expected. However, in SrCu02 the increase is 
very weak, which suggests an additional scattering mech- 
anism for phonons in this materialJ^ 



D. Theoretical treatment of the spinon-phonon 
scattering 

Expression [2] which describes the spinon-phonon scat- 
tering in our data has been empirically derived and is 
consistent with general considerations for Umklapp scat- 
tering. In the following we go beyond this empirical 
treatment and derive Zgp (T) from a spinon-phonon scat- 
tering theory presented in Ref. |4^ (labelled Ibp{T)). In 
this memory function approach it is possible to semi- 
analytically evaluate the temperature dependence of the 
mean-free path within the XY limit of the Heisenberg 
model assuming weak coupling. 

Applying a Jordan- Wigner transformation the XY 
Hamiltonian becomes, 

I I 

where — < — J/2 and A is the spin-phonon coupling 
constant in the spin-Peierls XY model. In this tight 
binding model the dispersion of fermions is given by 
e^, ~ — 2icos(fca) and the velocity at the Fermi wavevec- 
tor k ~ 7r/2a is equal to v — 2ta/h (a is the lattice con- 
stant). In the isotropic Heisenberg model the elementary 
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Figure 5: Fits to the mean free path according to the memory 
function approach explained in the text. The circles are the 
derived mean free paths of the spinon excitations in Sr2Cu03 
for 4N purity. The line represent a fit using Zap. The param- 
eters of the fit are cjq ~ 540 K and Iq = 1266 lattice spacings 
{lo = 4955 A). 



Figure 7: Comparison of the expressions for l^p only. The 
hexagonal shaped symbols represent an estimate for Isp de- 
rived from the magnetic mean free path Imag according to 
(T) = Z^ig (T) - IgK An lo = 1.56 ± 0.16 jim as approx- 
imated by the phenomenological model has been used. The 
line represents calculations of Isp with the memory function 
approach using a phonon frequency ujq = 540 K. 
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Figure 6: Fits to the mean free path according to the mem- 
ory function approach explained in the text. The pentagonal 
shaped symbols represent the derived mean free paths of the 
spinon excitations in SrCu02 for 4N purity. The line repre- 
sents a fit using an Up as defined by equation [4] The param- 
eters of the fit are luq — 540 K and lo = 3831 lattice spacings 
{lo = 15000 A). 



magnetic excitations - the spinons - have a velocity equal 
to Vsp = f To take into account the difference in ve- 
locities of the excitations between the XY and isotropic 
Heisenberg model we evaluate the mean free path using 
an effective hopping matrix element — tcff = f 7 • 

The mean free path is given by Igp (T) ~ VspT (T) where 
T is the characteristic scattering time. The phase space of 
the spinon-phonon scattering matrix elements determin- 
ing r is restricted by the energy-momentum conservation 
laws^ as, 

+ 7lq5{huj + Ek - €k+q + tlhJq)\u,^Q- (4) 

with the first term representing a phonon emission and 
the second, phonon absorption, fk = 1/(1 + e'^^''), 
Uq = \/[e^^^'i — 1) are the fermion and boson occupation 
factors respectively (/3 = l/fc^T")- 

It is clear from this formulation that low energy acous- 
tic phonons ujq ^ cq (c the sound velocity) cannot con- 
tribute to the scattering because the energy-momentum 
conservation laws cannot be simultaneously satisfied. We 
thus consider scattering by optical phonons of frequency 
Wo, as the effect of scattering by zone boundary acous- 
tic phonons (Umklapp scattering) is similar; we take a 
typical J ^ 2400 K. In the fits shown in Fig. [5] and 
Fig. [S] in the temperature range 30 K < T < 300 K we 
assume Z~^g (T) = Iq^ + (T) and optimize with re- 
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spect to Iq (the impurity scattering length) and ujq. As 
we do not know A, we normaUze the theoretical curve 
with respect to the experimental one at T = 300 K. For 
Sr2Cu03 (SrCu02) a good fit for the whole temperature 
range is found with = 1266 (^o = 3831) chain unit 
cells and ljq = 540 K for both compounds. A slightly 
better agreement at high temperatures can be achieved 
with a larger wq, although this leads to a strong devi- 
ation at low temperatures. The discrepancy at around 
100 K is due to the fact that the model does not accu- 
rately reproduce the slope of the experimentally deter- 
mined ^sp values below this temperature. This is illus- 
trated in Fig. [T] There, an l^p approximated from the 
measured data is compared to the Igp as obtained from 
the memory function approach. The approximation of 



^mig (T) - ^0 ^ where an l„ 



Isp is done by 1^^^ (T) ^ 

from equation [1] and the estimated Iq = 1.56 ± 0.16 pm 
by the phenomenological model is used. While the agree- 
ment between Igp and Igp is good at high temperatures, 
it gets poorer towards low temperatures. Especially in 
the double-logarithmic plot it can be seen that below 
100 K the model does not reproduce the slope of the ex- 
perimental estimate. It is interesting to note, that the 
fitted luq is somewhat lower but still of the same order of 
magnitude than typical frequencies of the optical Cu-0 
stretching mode^^— . In contrast, the phenomenologi- 
cal Isp according to equation [5] hints at a scattering by 
acoustical phonons as can be seen by the T* = 210 K. 

We mention that the exponential decay of the heat 
conductivity parallel to the spin chains Hb w Kmag ex 
exp(r*/T) is consistent with a theoretical treatment by 
Shimshoni et alj^i^. However, we do not observe Kc = 
Kph oc exp(2r^/T) as is expected in the same model. 



VI. SUMMARY 

We have investigated the spinon thermal conductiv- 
ity Kmag of high-purity single crystals of the single-chain 
5 = 1/2 AFM Heisenberg chain compound Sr2Cu03. We 
find that Kmag is strongly enhanced as compared to pre- 
vious results obtained on lower purity crystals. The anal- 
ysis of the data yields a very large low-temperature mean 
free path of ~ 0.5pm, corresponding to 1266 chain unit 
cells. Upon increasing the temperature towards room 
temperature, the mean free path decreases strongly and 
approaches that observed in lower purity samples. By us- 
ing a kinetic model we can decompose the mean free path 
into a temperature-independent spinon-defect scattering 
length Iq and a temperature dependent spinon-phonon 
scattering length l^p ^ ^ cxpT*/T with a characteristic 
energy scale ksT* for umklapp processes, and for low-T 
upon increasing the purity. 

By comparing the temperature dependence of the 
mean free path of Sr2Cu03 with that of SrCu02, we 
could show that the spin-phonon interaction, as ex- 
pressed by Igp is practically the same in both systems. 
The comparison of the empirically derived ^sp with model 





B in 10-=^^ K-^s^ 


A in 10"^^ 


L in 10"*m 


b 


2N 


2.25 


4.50 


7.50 


3.60 


4N 


3.45 


1.53 


2.33 


3.49 



Table L Fit parameters of a Callaway fit to Kc of Sr2Cu03, 
shown in the left panel of Fig. [T] with a — 3 and /? = 1. 



calculations for the spin-phonon interaction of the S — 
1 /2 AFM XY chain model yields a reasonable agreement. 
This agreement is very encouraging for further studies as 
an analysis of the full Heisenberg model might improve 
this agreement even more. 
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Appendix 

In order to model the phononic thermal conductivity 
perpendicular to the chain, a phenomenological model, 
devised by Callaway^, can be used. Although this model 
has undergone several revisions and extensions^"—, the 
main approach is to model Kph within the Debye approx- 
imation as^ 



2n'^Vph 



1) 



Tcdx. (5) 



Here x ~ fiw/kBT, lu is the phonon angular frequency, 
Q]j is the Debye temperature and Wph is the phonon ve- 
locity. In the Debye approximation the sound velocity is 
given as 



enk 



dkb 



(6) 



with N representing the number of elementary cells per 
unit volume, is a combined scattering rate, which is 
assumed to be the sum of all individual scattering rates 

r-' =Ts' +T^' +T^\ (7) 

where tb denotes boundary scattering, td point defect 
scattering and tjj Umklapp scattering. This separation 
is possible, as long as these scattering processes are inde- 
pendent of each other, which is the gist of Matthiessen's 
rule. In the context of the model, this is considered to 
be fulfilled since in different temperature regions different 
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scattering mechanisms are dominant. A Tc containing ex- 
pressions for all scattering processes can then be written 
as 

r-i = ^+^^4 + Sc.^rexp(-^), (8) 

with fit parameters L, A, B, b. The values of these pa- 
rameters are given in table |T] for the fit in the inset of 
Fig. [H The largest changes, when comparing the two 
purities, are for parameter A, which describes the concen- 



tration of point defects, and for parameter L, which de- 
scribes the boundary scattering. The difference in bound- 
ary scattering only indicates a difference in the sample 
geometries between the 2N and 4N samples. The de- 
crease of scattering by point defects upon an increase of 
purity however underpins the overall reduction of defects. 
It should be noted that in Ref. H two additional scatter- 
ing processes had to be used to describe the data. In 
our analysis which is focused on temperatures above the 
maximum of k these are not necessary. 
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